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Abstract 

Given a set I? of nonnegative integers, we derive the asymptotic number of graphs with a given 
number of vertices, edges, and such that the degree of every vertex is in T>. This generalizes existing 
results, such as the enumeration of graphs with a given minimum degree, and establishes new ones, such 
as the enumeration of Euler graphs, i.e. where all vertices have an even degree. Those results are derived 
using analytic combinatorics. 


1 Introduction 

1.1 Related works 

The asymptotics of several families of simple graphs with degree constraints have been derived. Regular 
graphs, where all vertices have the same degree, have been enumerated by Bender and Canfield (1978), 
graphs with minimum degree at least 5 by Pittel and Wormald (2003). An Euler graph, or even graph, is a 
graph where all vertices have an even degree. An exact formula for the number of such graphs, for a given 
number of vertices and without consideration of the number of edges, has been derived by Robinson (1969) 
and Mallows and Sloane (1975). In the present work, we generalize those results and derive the asymptotic 
number of graphs with degrees in any given set. 

A similar problem has been addressed with probabilistic tools by the configuration model, introduced 
independently by Bollobas (1980) and Wormald (1978). This model inputs a distribution F on the degrees, 
and outputs a random multigraph where the degree of each vertex follows F. The main difference with the 
model analyzed in this article is that the number of edges in the configuration model is a random variable. 
The link between both models is discussed in Section 4.1. For more information on the configuration model, 
we recommend the book of van der Hofstad (2014). 

Other related problems include the enumeration of graphs with a given degree sequence (Bender and 
Canheld (1978)), the enumeration of symmetric matrices with nonnegative coefficients and constant row 
sum (Chyzak et al. (2005)), and the enumeration of graphs with degree parities, investigated by Read and 
Robinson (1982). 

‘This work was partially founded by the Austrian Science Fund (FWF) grant F5004, the Amadeus program and the PEPS 
HYDrATA. 

tThis work was supported by Spains Ministerio de Ciencia e Innovacin under the project “Combinatoria, Teora de Grafos y 
Geometra Discreta” (ref. MTM2011-24097) 
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1.2 Model and notations 

A multiset is an unordered collection of objects, where repetitions are allowed. Sets are then multisets 
without repetitions. A sequence is an ordered multiset. We use the parenthesis notation (iti,...,u„) for 
sequences, and the brace notation {ui,... ,u„} for sets and multisets. Open real intervals are denoted by 
open square brackets ]a, h[. 

A simple graph G is a set V{G) of labelled vertices and a set E{G) of edges, where each edge is an 
unordered pair of distinct vertices. In a multigraph, the edges form a multiset and the vertices in an edge 
need not be distinct. An edge {u,w} is a loop ii v = w, a multiple edge if it has at least two occurrences in 
the multiset of edges, and a simple edge otherwise. Thus, the simple graphs are the multigraphs that contain 
neither loops nor multiple edges, i.e. that contain only simple edges. The set of multigraphs with n vertices 
and m edges is denoted by MG„_m, and the subset of simple graphs by SG„^rn- 

The degree of a vertex is defined as its number of occurrences in E{G). In particular, a loop increases its 
degree by 2. The set of multigraphs from MG„_m where each vertex has its degree in a set V is denoted by 
MG^^^i- The subset of simple graphs is SG^^Ji- The set V may be finite or infinite. We denote its generating 
function by 

Setc(x) = 

d^V 

For any natural number i, V — i denotes the set {d — i G Z>o | d € V}. In particular, observe that 
Set^(a;) = Setx)-i(a:). We also define the valuation r = min(I?) and periodicity p = gcdjfii — ^2 I di,d 2 G V} 
of the set (by convention, the periodicity is infinite when \'D\ = 1). 


2 Main Theorem and applications 


Our main result is an asymptotic expression for the number of graphs in SG^^^u when the number m of 
edges grows linearly with the number n of vertices. 

Theorem 1. Assume T) contains at least two integers, has valuation r = minjii G V} and periodicity p = 
gcdjdi — ^2 I di, d 2 G T)}. Let m, n denote two integers tending to infinity, such that 2m/n stays in a fixed 
compact interval o/]r, max(I?)[ and p divides 2m — rn, then the number of simple graphs in SG^^^i 


(2to)! Set-D(C)" P 

2™m! ^2TTnCfi'iO ™ 1 i ))> 

where (j){x) = C the unique positive solution of (j){C) = and W^(x) = If P 

does not divide 2m — rn, if2m/n < r or if2m/n > max(I?), then SG^^^i ts empty. 

When V = 1i>o, the degrees are not constrained, so SG^^^i = SG„_m- Using Stirling formula, it can indeed 

be checked that the total number of simple graphs with n vertices and m edges, is asymptotically 

equal to the result of Theorem 1 


.,2m 


( 2 to )! 


2™m! (2m)2"^e“^"‘-\/27r2m 




(l + 0(n-i)). 


Pittel and Wormald (2003) have derived the asymptotics of simple graphs with minimum degree at least 6. 
They used probabilitic and analytic elementary tools, in a sophisticated way. In the present paper, we have 
addressed the enumeration of a broader family of graphs with degree constraints, using more powerful tools 
(analytic combinatorics). For graphs with minimum degree at least S, the asymptotics derived in Theorem 1, 
for T) — matches their result. 

Euler graphs are simple graphs where each vertex has an even degree. An exact, but complicated, formula 
for the number of such graphs, for given number of vertices and without consideration of the number of edges. 
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has been derived by Robinson (1969) and Mallows and Sloane (1975). Applying Theorem 1, we are now able 
to derive the asymptotic number of Euler graphs with n vertices and m edges, when 2m/n stays in a fixed 
compact interval of R>o 


where ^(a;) 


(even) 


(2m)! cosh(C)" 2 
2-m! ^2TrnC(l>’{cf 


V 4m J 4m ^ 


(l + 0(n 1)), 


a;tanh(x) and tanh(C) = 2m/n. 


3 Proof of the result 

In this section we provide a proof for Theorem 1. The proof of all lemmas and theorems are moved to the 
appendix. 


3.1 Preliminaries 

3.1.1 Multigraph model 

The main model of random multigraphs with n vertices and m edges is the multigraph process, analyzed 
by Flajolet et al. (1989) and Janson et al. (1993). It samples uniformly and independently 2m vertices 
{vi,V 2 , ■ ■ ■, V 2 m) in {1) • ■ •, n}, and outputs a multigraph with set of vertices {1,..., n} and set of edges 

{{v 2 i-i,V 2 i} I 1 < i < m}. 

Given a simple or multi graph, one can order the set of edges and the vertices in each edge. The 
result is a sequence of ordered pairs of vertices, that we call an ordering of G. Let orderings(G) denote the 
number of such orderings. For example, the multigraph on 2 vertices with edges {{1,1}, {1, 2}, {1, 2}} has 12 
orderings, amongst them ((1, 2), (1,1), (2,1)). For simple graphs, the number of orderings is equal to 2”*m!, 
because each edge has two possible orientations and all edges can be permuted. For non-simple multigraphs, 
orderings is smaller. Flajolet et al. (1989) and Janson et al. (1993) introduced the compensation factor k{G) 
of a multigraph G with m edges, defined as 

_ orderings(G) 

“ -dlJ—i-■ 

2”^m! 

The compensation factor of a multigraph is 1 if and only if it is simple. 

Observe that in the random distribution induced by the multigraph process, each multigraph receives a 
probability proportional to its compensation factor. Therefore, when the output of the multigraph process is 
constrained to be a simple graph, the sampling becomes uniform on The total weight of a family iF 

of multigraphs is the sum of their compensation factors. For example, the total weight of MG„_m is equal 

2m 

to When F contains only simple graphs, its total weight is equal to its cardinality. 


3.1.2 Analytic tools 


Our tool for the analysis of graphs with degree constraints is analytic combinatorics, as presented by Flajolet 
and Sedgewick (2009). Its principle is to associate to the combinatorial family studied its generating function. 
The asymptotics of the family is then linked to the analytic behavior of this function. 

In the analysis of a graphs family F with analytic combinatorics, the main difficulty is the fast growth 
of its cardinality, which often implies a zero radius of convergence for the corresponding generating function 




z\V{G)\ 


This feature drastically reduces the number of tools from complex analysis that can be applied. Graphs with 
degree constraints are no exception, but our approach completely avoid this classic issue. In fact, the only 
analytic tool we use is the following lemma, a variant of (Flajolet and Sedgewick, 2009, Theorem VIII.8). 
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Lemma 2. Consider a non-monomial series B{z) with nonnegative coefficients, analytic on C, with valuation 
r = min{n | [z"]i?(z) 7 ^ 0} and periodicity p = gcd{n | 7 ^ 0}. Let (j){z) denote the function 

and K a compact interval of the open interval ]r, liina;_>.oo Let N, n denote two integers tending to 

infinity while N/n stays in K, and let f denote the unique positive solution of 4>{C,) = N/n. Finally, consider 
a compact Y and a function A{y,z), on Y x C, such that for all y in Y, the function z 1 —> A{y,z) is 
analytic on C and A{y, (ffi) is nonzero. Then we have, uniformly for N/n in K and y in Y, 


[z^\A{y,znB{z) 


pMvX B(c) lU if p divides N — nr, 

■y/2imC<^'(C) 

0 otherwise. 


3.2 Multigraphs with degree constraints 

The work of Flajolet et al. (1989) and Janson et al. (1993) demonstrates that multigraphs are more suitable 
to the analytic combinatorics approach than simple graphs. Moreover, the results on multigraphs can usually 
be extended to simple graphs. Following this observation, multigraphs are analyzed in this section, before 
turning so simple graphs in Section 3.3. 


3.2.1 Exact and asymptotic enumeration 

We derive an exact expression for the number of multigraphs with degree constraints in Theorem 3, then 
translates it into an asymptotics in Theorem 4. 

Theorem 3. The total weight of all multigraphs in is 


The proof of this theorem is elementary by the definition of the compensation factor. Now applying 
Lemma 2 to the exact expression, we derive the asymptotics of multigraphs with degree constraints. Let us 
first eliminate three simple cases. 

• When V contains only one integer V = {d}, is the set of d-regular multigraphs. The total 

weight of is then 0 if 2 m 7 ^ nd, and otherwise. 

• The sum of the degrees of the vertices is equal to 2m, so i® empty when 2m/n < min(I?) or 

2m jn > max(I?). 

• The periodicity p of I? is equal to gcdjd— r | d G V}. For each vertex w of a multigraph from MG^^^d 

it follows that p divides deg(u) — r. By summation over all vertices, we conclude that if p does not 
divide 2 m — nr, then the set i® empty. 

The two last points obviously hold for SG)j®2i- 

Theorem 4. Consider a set B C Z>o of size at least 2. Let r = min(I?) denote its valuation and p = 
gcdjdi — d 2 I di,d 2 G B} its periodicity. Let m, n denote two integers tending to infinity, such that 2m/n 
stays in a fixed compact interval of the open interval ]r,m.ayi(B)[, andp divides 2m —rn, then the total weight 
o/MG^®^ is equal to 


E = 

ggmg).®,). 


(2m)! p Setx)(C)” 
2’^m! 


(l + 0(n-i)) 


where 4>{x) = ^ ^ unique positive solution of (/{(,) = 


4 











2 1 3 



Figure 1: Four examples of multigraphs from 


3.2.2 Typical multigraphs with degree constraints 

Let us recall that an edge is simple if it is neither a loop nor a multiple edge. Before turning to the 
enumeration of simple graphs with degree constraints, we first describe the behavior of non-simple edges in 
a typical multigraph from proofs are given here, as stronger results will be derived later. 

Using random sampling, we observe that in most of the multigraphs from all non-simple edges 

have low multiplicity and are well separated. This motivates the following definition. A multigraph 
from is in if all its non-simple edges are loops or double edges, and each vertex belongs to 

at most one loop or (exclusive) one double edge. Let \E\e denote the number of occurrences of the element e 
in the multiset E. Formally, is characterized as the set of multigraphs G from MG^^^ such that 

for all vertices u, v, w, we have 

< 1, = |A(G)|{,,^} =2 ^ u = w, 

< 2, {v,v} e E{G) Ww, < 1. 

The complementary set, MG^®^ \ is denoted by and illustrated in Figure 1. 

3.3 Simple graphs with degree constraints 

We introduce the notation SG^^^i the set of simple graphs with n vertices, m edges and all degrees in T), 
i.e. multigraphs from MG^^ that contain neither loops nor multiple edges. The enumeration of simple 
graphs with degree constraints is derived in Theorem 1. First, in Section 3.3.2, we describe an inclusion- 
exclusion process that outputs | SG^^^i I when applied to In Section 3.3.3, this process is then 

applied to MG^^, and the error introduced is proven to be negligible in Section 3.3.4. 

In order to forbid loops and multiple edges in multigraphs from MG^^, we introduce the notion of 
marked multigraphs. 

3.3.1 Marked multigraphs 

A marked multigraph G is a triplet (U(G), i5(G), A(G)), where U(G) denotes the set of vertices, E{G) the 
multiset of normal edges, and E{G) the multiset of marked edges, where both normal and marked edges 
are unordered pairs of vertices. We say that a marked multigraph G belongs to a family E of (unmarked) 
multigraphs if the unmarked multigraph {V{G),E{G) U E{G)) is in E. 

We now extend to marked multigraphs the definitions of degree, orderings and compensation factors, 
introduced for multigraphs in Section 3.1. The degree of a vertex from a marked multigraph G is equal 
to its number of occurrences in the multiset E{G) U E{G). An ordering of a marked multigraph G with 
TO = |£'(G)| -I- |A(G)| edges is a sequence 

S = {{vi,Wi,ti), ..., {Vm,WmEm)) 

from (U(G) x U(G) x {0,1})™ such that the multiset {{vi,Wi} \ {vi,Wi,0) G S'} is equal to E{G), and 
the multiset {{vi,Wi} \ {vi,Wi,l) G S} is equal to E(G). The number of orderings of a given marked 
multigraph G is denoted by orderings(G), and its compensation factor is 

_ orderings(G) 

Klijr) — ~ ~ ^ , 
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For example, consider the marked multigraph G with 


ViG) = {1,2}, E{G) = {{1,2}}, EiG) = {{1,2}, {1,2}}. 

Its number of orderings is 24, and therefore its compensation factor is k{G) = 1/2, whereas it is 1/6 for G 
without the marks, 

V{G) = {1,2}, EiG) = {{1,2},{1,2},{1,2}}. 

In the following, we will consider families of marked multigraphs where the marked edges are loops or 
multiple edges. Given a marked multigraph G, then £{G) denotes the number of loops in E{G), and k{G) 
the number of distinct edges from EiG) that are not loops. The generating function of a family E or marked 
multigraphs is 

Fiu,v) = KiG)u^^^\^^^\ 

GeJ^ 


3.3.2 Inclusion-exclusion process 

In this section, we build an operator Marked that inputs a family of multigraphs and outputs a family of 
marked multigraphs. It is designed so that the asymptotics of its generating function Markodj^^co) (m, v) 

is linked to the asymptotics of | SG^^^i I- order to justify the construction, we first introduce the opera¬ 
tors Marked^^^ and Marked^^^ 

First marking. If we could mark all loops and multiple edges from MG^®2i) the enumeration of simple 

graphs with degree constraints would be easy. Indeed, given a family E of multigraphs, let Marked^^ denote 
the marked multigraphs from E with all loops and multiple edges marked. Since the simple graphs are the 
multigraphs that have neither loops nor multiple edges, we have 

Marked^|^j^^(0,0)= ^ ^ ^ 

g^mgEX g^sgEX 

which is equal to | SG^^^i 1; because simple graphs have a compensation factor equal to I. Unfortunately, we 
do not have a description of this family in the symbolic method formalism. 

Second marking. The inclusion-exclusion principle advises us to mark some of the non-simple edges. 
Let Marked^"^ denote the set of marked multigraphs G from E such that each edge from EiG) is either a 
loop, or has multiplicity at least 2 in EiG) and does not belong to EiG). This construction implies the 
relation 

Marked^^( m, v) = Marked^^('u -I- l,u -I- 1), 

and therefore 

|SGl®j,|=Marked|^|^<^^(-l,-I). 

( 2 ) 

The natural idea to build a marked multigraph G from Marked^^(,p) is to first choose some loops and 

multiple edges to put in EiG), then complete EiG) with unmarked edges, which may well form other loops 

( 2 ) 

and multiple edges, in a way that ensures G € Marked^^(j 3 ) . However, the description of the set of marked 
edges is complicated, because of the numerous possible intersection patterns. 
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Third marking. We have seen in Section 3.2.2 that in most of the multigraphs from MG^^, non-simple 
edges do not intersect. This motivates the following definition. Given a set T of multigraphs, let Marked(J^) 
denote the set of marked multigraphs from T such that each vertex is in exactly one of the following cases: 

• the vertex belongs to no marked edge, 

• the vertex belongs to one marked loop and no other marked edge. 


• the vertex belongs to two identical marked edges and no other marked edge. 

Therefore, each marked edge is a loop of multiplicity 1 or a double edge. This marking process links the 
multigraphs from MG^^\ defined in Section 3.2.2, to the simple graphs with degree constraints. 

Lemma 5. The value Markodj^^^cu,*) (— 1, — 1) is equal to the number of simple graphs in SG|j^2i- 

Applying the operator Marked to the decomposition 

^n,m ^n,m ^ ^n,m ’ 


we find 


Markedj^Q{c) (u, v) = Markedj^^cD,.) {u, v) + Markedj^Q{®,o) (u, v) 
which implies, after evaluation at {u,v) = (- 1 ,- 1 ) and reordering of the terms, 

I I = ^^arked^Qjin (-1, -1) - Markedj^^ji^i) (-1, -1). 

We compute the asymptotics of Markodj^^cc) (—1, —1) in Section 3.3.3, and prove that Markedj^^cc.o) (— 
is negligible in Section 3.3.4. 


1 ,- 1 ) 


3.3.3 Application of the inclusion-exclusion process to all multigraphs with degree constraints 

We start with an exact expression of Markedj^^cx,) (u, v) in Lemma 6 , then derive its asymptotics in Lemma 8 . 
Lemma 6. We have the formal equality 


MarkedMCjW (u, ?;) = 


,m,2fc+£ 

. fc,£>0 


{uW2r{xYY {vWir{x)Y 


k\ 


i\ 


Set-p(a:)^ 


where an,m,j = 0 when j is greater than min(n, m), otherwise 

n\ ml {2m — 2j)\{2mY^ 

{n — j)\ni {m — j)\mt ( 2 m)! ’ 

’ Am Seti,(x) ■ 

The proof is constructive by considering all the disjoint sets of vertices where we can put a loop or a 
double edges. We observe that when 2k + i is fixed while n,m tends to infinity, then an,m, 2 k+i tends to 1. 
The double sum can then be approximated by an exponential, and it is tempting to conclude 

Markedj^Q(i3) {u,v) ^ Setx>(a:)". 

The next lemma formalize this intuition. A multivariate generating function f{xi ,..., x„) is said to dominate 
coefficient-wise another series g{xi ,..., Xn) if for all fci,..., fcn > 0 , 


[xY ■ • • xY-]g{xi, ...,x„) < [xY ■ ■ ■ xY-jfixi ,..., Xn). 
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Lemma 7. When mjn stays in a fixed compact interval o/]R>o, there is an entire bivariate analytic func¬ 
tion C{u,v) such that, for n large enough, ^C{u,v) dominates coefficient-wise 


eu+v 


E 

k,£>0 




h P 
U V 


( 1 ) 


We can now derive the asymptotics of Markadj^^cc) {u,v). As observed in the discussion preceding 
Theorem 4, the result is trivial when V contains only one integer, when 2mjn is outside [min(I?), max(I?)] 
and when p does not divide 2m — min(I?)ri. 

Lemma 8. Assume TA has size at least 2, valuation r and periodicity p. Let m, n denote two integers tending 
to infinity, such that 2m jn stays in a fixed compact interval of]r, max(I?)[ and p divides 2m — rn. When u, 
V stay in a fixed compact, then 


Markedj^^CD^ {u, v) 


(2m)! Setx)(C)” p 

2 ™m! ^27rnC((>'(^ 




(l + 0(u-i)), 


( 2 ) 


where IT^ (cc) 


n a;^Setp_ 2 (x) j/ _ x Setp-i(a:) 
4m Setx)(ai) ’ ' Setx>(ic) 


and (j){C,) 


2m 

n 


The proof is a consequence of Lemma 2, Lemma 6 and Lemma 7. 


3.3.4 Negligible marked multigraphs 

Recall that denotes the set \ In Lemma 10, we prove that Markedj^Q(i3,o) (-1,-1) 

is negligible. To do so, we first bound Markedi^(l, 1) for a family R of marked multigraphs from 
with mandatory edges. 

Lemma 9. Let ei,..., denote j edges on the set of vertices {1, ..., n}, and R the set of multigraphs from 
MG^^2 i contain those edges, with multiplicities (i. e. an edge with k occurrences in the list has at least k 
occurrences in the multiset of edges of the multigraph) 


R = 


|g € VI < i < j, Ci S E{G) with multiplicities'^. 


Assume R contains at least two integers and has valuation r. Let m, n denote two integers tending to 
infinity, such that 2m/n stays in a fixed compact interval o/]r, max(I?)[, then 

Markedi^(l, 1) = 0{n~^ Markedj^^cc) (1,1)). 


Figure 1 displays four multigraphs from Actually, any multigraph from contains one 

of those four graphs as a subgraph, and this property can be described in terms of mandatory edges. In the 
following lemma, we use this fact to bound Markedj^gcc.o) (-1,-1)- 

Lemma 10. Assume R contains at least two integers, has valuation r and periodicity p. Let m, n denote 
two integers tending to infinity, such that 2m jn stays in a fixed compact interval o/]r, max(I?)[, and p 
divides 2m — nr, then 


Markedj^Q(i3,o) (—1, —1) = O (ji ^ Markedj^gcc) (—1, —1)^ . 

The intuition supporting this proof is that a multigraph G belongs to if and only if it contains a 

vertex v that is in one of the four configurations depicted in Figure 1. According to Lemma 9, multigraphs 
from MG^®^ that contain those subgraphs have a negligible total weight. Now we have all the ingredients 
to prove Theorem 1. 







Proof of Theorem 1. In Lemma 5, we have proven that the number of simple graphs in is equal to 

Markedj^Q(i7,.) (—1, —1). By a set manipulation, this quantity can be rewritten 

Markedj^QOn (-1,-1) - Markedj^Q^c^o) (-1,-1), 

where \ Replacing the second term with the result of Lemma 10, we obtain 

I SGi^i, I = Markedj^g^^(-l,-l)(l + Oin-^)). 

Finally, the asymptotics of Markedj^Q(c) (-1,-1) has been derived in Lemma 8. □ 

4 Random generation 

In order to keep a combinatorial interpretation, we focused on generating functions Setx) {x) with coefficients 
in {0,1}. Our results hold more generally for any generating function D{x) with nonnegative coefficients and 
large enough radius of convergence (so that the saddle-point from Lemma 2 is well defined). Multigraphs 
are then counted with a weight that depends of the degrees of their vertices 

weight(G) = k{G) deg('(;)![a;‘^®®^’'^]D(a:). 

vev(G) 

The present work has been guided by experiments on large random graphs with degree constraints. We 
used exact and Boltzmann sampling (Duchon et al. (2004)). Observe that to build a random simple graph 
from SGi^i,, one can sample multigraphs from reject until the multigraph is simple. As a 

consequence of Theorem 1, the expected number of rejections is e (using the notations of 

the theorem). 


4.1 Boltzmann sampling 


The construction of the Boltzmann algorithm is straightforward from Theorem 3. To build a random 
multigraph with degrees in V, n vertices and approximately m edges, the algorithm first computes a positive 
value X, according to the number of edges targeted. It then draws independently n integers (di,..., d„), 
following the law 


P(d) 


{[z‘^]D{z)) x‘^ 
D{x) 


(3) 


with D{x) = Setx)(a:). If their sum is odd, a new sequence is drawn. Otherwise, the algorithm outputs a 
random multigraph with sequence of degrees (di,..., d„). To do so, as in the conhguration model (Bollobas 
(1980), Wormald (1978)), each vertex Vi receives di half-edges, and a random pairing on the half-edges is 
drawn uniformly. 

Therefore, the random distribution induced on multigraphs by the Boltzmann sampling algorithm is iden¬ 
tical to the conhguration model. Conversely, given a probability distribution on Z>o, one can choose D{x) 
so that the distribution is equal to the one described by Equation (3). Thus, we expect random multigraphs 
from the conhguration model and multigraphs with degree constraints to share many statistical properties. 


4.2 Recursive method 

For the sampling of a multigraph in the generator hrst draws a sequence of degrees, and then 

performs a random pairing of half-edges, as in conhguration model and the Boltzmann sampler. Each 
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sequence (di,... ,d„) from is drawn with weight n«=i In th® first step, we use dynamic pro¬ 

gramming to precompute the values {Sij)o<i<n,o<j< 2 m, sums of the weights of all the sequences of i degrees 
that sum to j 


= 




n 


di,...,diGT> v—1 
diH- ^di=j 


dj’ 


using the initial conditions and the recursive expression 


[1 if (z,j) = (0,0), 

“ 1 ^ if * = 0 and j ^ 0, or if j < 0, 

[j2dev otherwise. 


After this precomputation, we generate the sequence of degrees as follows: first we sample the last degree d„ 
of the sequence according to the distribution 


P(d„ = d) 


Sn—l,2m—d 

d^-Sji^2m 


then we recursively generate the remaining sequence (di,..., d„_i), which must sum to 2m — dn- Once the 
sequence of degrees is computed, we generate a random pairing on the corresponding half-edges. 


5 Forthcoming research 

The results presented can be extended in several ways. The case where 2m jn tends to min(I?) or max(I?) 
could be considered. For example, Pittel and Wormald (2003) have derived, using elementary tools, the 
asymptotics of graphs with a lower bound on the minimum degree when m = 0(n log(n)). This extension 
would only require to adjust the saddle-point method from Lemma 2. 

We have also derived results on the enumeration of graphs where the degree sets vary with the vertices. 
The model inputs an infinite sequence of sets (I?i, ■ • ■) and output graphs where each vertex v has 

its degree in Vy. The techniques presented in this paper can be extended to this case, if some technical 
conditions are satisfied, such as the convergence of the series n~^ log(Set-D„ (x)). This extension will be 
part of a longer version of the paper. Two examples of such families are graphs with degree parities (Read 
and Robinson (1982)), and graphs with a given degree sequence (Bender and Canfield (1978)). 

We believe that complete asymptotic expansion can be derived for graphs with degree constraints. This 
would require to apply a more general version of Lemma 2, such as presented in Chapter 4 by Pemantle and 
Wilson (2013), and we would have to consider more complex families than 

The asymptotics of connected graphs from when m — n tends to infinity has first been derived by 

Bender et al. (1990). Since then, two new proofs were given, one by Pittel and Wormald (2005), the other 
by van der Hofstad and Spencer (2006). The proof of Pittel and Wormald relies on a link between connected 
graphs and graphs from a particular family of graphs with degree constraints (graphs with degrees at least 2). 
In de Panafieu (2014), following the same approach, but using analytic combinatorics, we obtained a short 
proof for the asymptotics of connected multigraphs from MG„^m when m — n tends to infinity. We now plan 
to extend this result to simple graphs, and to derive a complete asymptotic expansion. 

In this paper, we have focused on the enumeration of graphs with degree constraints. We can now start 
the investigation on the typical structure of random instances of such graphs. An application would be the 
enumeration of Eulerian graphs, i.e. connected Euler graphs. 

Finally, the inclusion-exclusion technique we used to remove loops and double edges can be extended to 
forbid any family of subgraphs. 
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A Proofs 

In this appendix, we include the proofs of the lemmas and theorems. 


Proof of Theorem 3. By definition of the compensation factor, the number of multigraphs of the theorem is 
equal to 

E orderings(G). 

G^MG'ff’X 


Let us consider an ordering 

((^^ 1 , ^' 2 ), {V2,V3), iv2m-l,V2m))- 

of a multigraph G from For all 1 < i < n, let Pi = {j \ vj = i} denote the set of positions of the 

vertex i in this ordering. Since the vertices have their degrees in V, each Pi has size in T). This implies a 
bijection between 

• the orderings of multigraphs in 


• the sequences of sets (Pi,..., P„), where the size of each set is in V, and (Pi,..., P„) is a partition 
of {1,..., 2m} {i.e. the sets are disjoints and Ur=i -P* = {1 j • • ■: 2m}). 

We now interpret (Pi,... ,P„) as a sequence of sets that contain labelled objects and apply the Symbolic 
Method (see Flajolet and Sedgewick (2009)). The exponential generating function of sets of size in V is 
Setx)(a;). The bijection then implies 


orderings(G) = (2m)![a:^™] Set-D(a;)", 

g&mgXX 


and the theorem follows, after division by 2™ml. □ 

Proof of Lemma 5. As explained in the paragraphs First markink and Second marking of Section 3.3.2, 
the following relations hold 


Marked;^J^,^,(0,0) = |SGi®2,|, 

Marked^^^ (it, v) = Marked^^^ (it + 1, u + 1). 

Furthermore, by construction of , we have 

Markedj^QCD^.) (it, v) = MarkedJ^^ji,,., (it, v), 


so Markedj^Q^c^^,) (-1, -1) = | SG^^^. I- ^ 

Proof of Lemma 6. To build an ordering of a multigraph from Markedj^Q(x>) with 2k vertices in marked 
double edges and £ vertices in marked loops, we perform the following steps: 

1. choose the labels of the 2k vertices that appear in the marked double edges, and the £ vertices that 
appear in the marked loops. There are ( 2 ^, ^ X- 2 k-^ such choices. 

2. choose the distinct k edges of distinct vertices, among the chosen 2k vertices, that will become the 
marked double edges. There are such choices. 

3. order the 2k marked double edges and the vertices in each of them. There are ways to order 

them. 
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4. order the i loops. There are ^\ ways to do so. 

5. choose among the m edges of the final ordering which ones receive marked loops and which ones receive 
marked double edges. There are i rT- 2 k-e) choices. 

6. to fill the rest of the final ordering, build an ordering of length 2m — Ak — 2£ where the 2k vertices 

that belong to marked double edges and the i vertices that appear in marked loops have degree 
in I? — 2, while the other n — 2k — I vertices have degree in T). The number of such orderings is 
(2m - 4fc - 2£)![a;2™-4fc-2^] Set75(a;)"-2'=-^. 

This bijective construction implies the following enumerative result 




GGMarked(MGC(l!,) 
“ 2" 


m 

,m — 2k — £ 


n \ (2fc)!(2fc)!4fe / 

'^m\ ^ \2k,£,n — 2k — £/ 2^k\ 2^ \2k,£,i 

(2m -4k- 2£)![a;2™-4'=-2^] 

After simplification, this last expression can be rewritten 


Markedjyji^CD^ (li, v) — ^ j [x ] | ^ ^ nn,m, 2 fc+<' 


(uW2L(xyf {vW:l{x)Y 


, fc.«>0 


k\ 


£\ 


Set-D(x)"'. 


□ 


Proof of Lemma 7. Developing the exponential as a double sum 


e“+" = 


E 

k,e>o 


h f 
U V 

kill' 


the result can be rewritten 

\4 — an,m,2k+(.\ , r k \ 

"- k\£\ -- ^ 

for all k, £. We prove that when n is large enough, we have 

^jl - a n,-m,2k+l\ ^ G I ^ \ (2fc + 

” ^ k\£\ - \ m) y/kW. 


for all fc,£ > 1. Since the right-hand side are the coefficients of a function analytic on C^, this will conclude 
the proof. 

Let bn,j denote the value 11^=0 ~ n) > observe that an,m,j is equal to bnjbm,j/b 2 m, 2 j- Since 

bn,j < 1, if {cn,j) denotes a sequence such that < bnj for all then CnjCmj < o,n,m,j < 

which implies 


|1 — an,Tn,2k+e\ , ^^'^i'^2m,4k+2e 


k\£l 


< n- 


I5 1 Cn,2k+^^m,2k+Y) 

~m\ 


(5) 


We now prove that Equation (4) holds both for 2k + £ < y/rnl2 and for 2k + £> y/rnf2. 

• 2 

Case 2k + £ < y/rnf2. We prove by recurrence on j that bnj > 4 — ^ - The recurrence is initialized with 
bn,o = 1- Assuming it is satisfied at rank j, then 


bn,j+l 





> 


O' + I)" 
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which concludes the proof of the recurrence. This implies, using Inequality (5) 


|1 Oin,m,2k-\-£\ ^ ^ ( f 

"- m -— 


{Ak+21Y 

2m 


- 1 , 1 - 1 - 


{2k + £Y 


1 - 


{2k + £f 


m 


Since 2k + £ < C™/2, the first argument of the maximum function is at most 1. The second argument is 
smaller than {n~^ + m~^){2k + £)^. Therefore, we have 


|1 ^n,m^2k+t\ ^ 

m\ ' - 


V 771/ 


n \ {2k + £Y 


k\£\ 


and Inequality (4) is satisfied. 

Case 2k + £ > •y/w/2. We first prove bnj > e~^. To do so, we apply a sum-integral comparison in the 
expression 


/-i 


i=0 


^og{bn,j) = T ^ - / log (l - 0 = -(?7 - j) log f 1 - ^ j - j > -j. 


Inequality (5) then implies 

^ ^ max - 1,1 - < 

fpl'C. /vil-t-. ' ' 


n 

VkW. VkW. ' 


We now prove that nj\/ld£\ is smaller than 1 for n large enough. Indeed, 2k+ £ > yTn/2 implies max(fc,£) > 
y/m/8, so 


< 


< 


\/kl£l Cmax(A;,^)! (v^/8)!’ 

and since m/n stays in a compact interval of IEl>o, this last term tends to 0 with n. We then conclude 


|1 ^n,m^2k+£\ ^ e 

n -r^TT-^- < 


4k+2i 


k\£\ - CM 

for n large enough, so Inequality (4) is satisfied. 

Proof of Lemma 8. We start with the expression of Markedj^^^cD) (u, v) derived in Lemma 6 


□ 


( 2??? 1 ^ 

Marked^ci-)(7r,.) = ^C^-] 


Using the notation 


nWjTL (xy+vWjn (x) 


^ {uW^{xf)>^ ivW^{x)y ^ 

/ ^ ^n,m,2k+£ ^ | | oetx)(^) 


. k,e>o 


k\ 


£\ 


,m,2/c+£ 

fc,£>0 


(wIUtl(x)^)^ {vWj^{x)Y 


k\ 


A{x) = e 

this implies 

(2m)! |. 2mi uWr!,(x)^+vWr!,(x) 

2™m!^ ^ 2™m!' 

Observe that Wil{x) has valuation 0 and period p. According to Lemma 2, we have 

(X) ^ Setp(C)" ^ (cf+vW^ (C) ^ 0{n-^)), 


£l 


(2w)!r^2mi 


Setx)(a;)" - Marked cm {u,v) = ——'.[x "^]A{x) Set-Dix) 


2"*m!' 


2”‘m! C"™ C2^nCC(C) 
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so the demonstration is complete if we prove 




,_i Seti,«)- 


2™m!^ J \ ^ ^ ^ 2"*m! 

The Taylor coefficients of Wil (x) need not be positive, so we introduce the entire function 


n>0 


which dominate Wjl{x) coefficient-wise. By application of Lemma 7, (x)^,xIT jl( x)) dominates 

coefficient-wise A{x), and therefore 


g^[x2“]^(a:)Set^(x)" 

2™m! 


- ^( 2 :)) Seti,(x) 

2"*to! n "* "* 


Finally, according to Lemma 2, we have 


[x^™] ^ C(wfF^ (x)2, (x)) Seti,(x) 

2"*to! n 


(2m)!^/ _i Setp(Cr \ 
2™m! V y 


□ 


Proof of Lemma 9. Let i? denote the set of multigraphs from with j distinguished mandatory edges 


ei = {vi,V2}, ■ ■ ■ ,ej = {v2j-i,V2j}. 


Given an ordering of a multigraph from R, we can distinguish the first occurrences of the mandatory edges, 
in order to obtain the ordering of a multigraph from R. Therefore, the number of orderings of multigraphs 
from R is at most equal to the number of orderings of multigraphs from R. Dividing by 2™m!, this implies 

^ n{G) < Y, «(G), 

Gefl Geii 


so Marked_R(l, 1) < Marked^(l, 1). 

Let W denote the fixed set of vertices that appear in the mandatory edges, and for all w G W, let 
denote the number of occurrences of the vertex w in the mandatory edges 

dw = \{i \Vi = ui}|. 

Let also Gn}m denote the set of multigraphs with n vertices and m edges, where each vertex w from the 
mandatory edges has degree in I? — and the other vertices have degrees in T). To construct an ordering 
from a multigraph in Marked^, we choose the j positions of the mandatory edges among the m positions 
available, the order of the vertices in those edges, and mark or not each of them. Then the rest of the 
ordering is filled with an ordering from Marked^cd) . Therefore, the number of orderings from Marked^ is 
at most 

- j)! Marked^(d) (1,1). 

Dividing by 2'"m! and using the fact that j is fixed, we obtain 

Marked^(l, 1) = O ^Markedg(d) (1,1)^ . (6) 

Following the steps of Lemma 6, Markedg,(d) (1,1) is smaller than or equal to 


(2m-2j)! 

2"*“J(m — j)\ 


X 


2m-2j-l 1 „ 

J I / (^n,m—j,2k-\-i 

\ k,i>0 


|^^(x)2/c+G 


Set-D_d„ 



Seti,(x)"-l^l. 


15 














An application of the same argument as in the proof of Lemma 8 leads to 


Marked^^.,(„„^_^.)(l,l) 


(2m-2j)! Setp(C)"-l^l \ 


Since |VL| and j are fixed, this implies, using Lemma 8 , 

( 2 m — 2 'iV 2 '"m' / \ 

Marked^^.)^_^ (l, 1) = 2m-j (Marked^g^m (1,1)) • 

Simplifying and injecting this relation from Equation ( 6 ), we obtain 

Marked^(l, 1) = 0{n~^ Markedj^^^cc) (1, !))■ 

□ 

Proof of Lemma 10. By definition, a multigraph G belongs to if and only if it contains a vertex v 

that is in one of the following configurations: 

1 . the loop {v,v} appears at least twice in E{G), 

2 . there is a vertex u such that the edge {u,v} appears at least three times, 

3. there is a vertex u such that {v,v} is in EifG) and {u,v} appears at least twice, 

4. there are vertices u and w such that {u,v} and {v,w} both appear at least twice. 

Let .Ri (resp. R 2 , R 3 , R 4 ) denote the set of multigraphs from that contain a vertex in configuration 1 

(resp. 2, 3, 4). We then have 

= Ri U R2 U R3 U R4. 

Let also Ri, R 2 , R 3 and R 4 denote four subsets of such that 

1 . the multigraphs from Ri contain two occurrences of the loop { 1 , 1 }, 

2 . the multigraphs from R 2 contain three occurrences of the edge { 1 , 2 }, 

3. the multigraphs from R 3 contain an occurrence of {1,1} and two occurrences of { 1 , 2 }, 

4. the multigraphs from R 4 contain two occurrences of {1,2} and two occurrences of {1,3} 

(see Figure 1). Given the symmetric roles of the vertices, the number of orderings from multigraphs in Ri 
(resp. R 2 , R 3 , R 4 ) is lesser than or equal to n times (resp. n^, n^, rt') the number of orderings from 
multigraphs in Ri (resp. R 2 , R 3 , R 4 ). This implies 


Marked^^ (1,1) < n Marked^jj (1,1), 
Marked^^(l, 1) < Markedij^ (1,1), 

Marked^^(l, 1) < Marked_R 3 (1,1), 

Marked^^(l, 1) < Marked_R„ (1,1), 


so 

Markedj^Q(i3,o) (1,1) < n Marked^jj (1,1) + Markedi^^ (1,1) + Marked^^g (1,1) + Markedi^3(l, 1). 

The multigraphs from Ri (resp. R 2 , R 3 , R 4 ) have 2 mandatory edges (resp. 3, 3, 4). Four applications of 
Lemma 9 lead to 

Markedj^g(T),o) (1,1) = 0{n~^) Markodj^^cc) (1,1). 
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Finally, according to Lemma 8, 

Markedj^Q^m (1,1) = O (^Markedj^^^m (-1, 
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